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SUMMARY 

. Analysis of the supersonic flow In two two-dimensional hi^- 
solldlty cascades and In a partly supersonic symmetrical nozzle shows 
that there Is, In general, significant deviation of the mean streamline 
shape from that of the mean hlade line and that the effect of hlade 
thickness and hlade curvature on the specific mass flow along the mean 
streamline Is to Increase the specific mass flow along the mean stream- 
line about 9 percent above that given by a one-dimensional estimate. 

In order to determine these effects more accurately for ttirbomachines 
of arbitrary hub and casing shapes to be used for the three-dimensional 
through- flow calculation, a method is developed for the determination of 
the supersonic flow along stream surfaces of revolution in turbomachines. 
In this method, the shapes of the stream surfaces are arbitrary, and the 
method also takes into account the distance between adjacent stream 
surfaces, which varies along the flow path. Thus, the method can be 
applied to t-urbomachines with arbitrary hub and casing shapes. 

In addition to their use for direct problems, these equations can 
be used to design blade elements in supersonic flow along an arbitrary 
stream filament of .revolution in turbomachines . 


INTRODUCTION 

Recent investigations of the applicability of supersonic flow in 
compressors have shown the desirability of using such flow to increase 
the pressure ratio per stage (references 1 and 2) . For axial-flow 
compressors having blades with short radial length, the over-all 
performance analysis is often based on a one- dimensional approximation, 
in which only average values in the channel are considered (for exan^le, 
references 3 and 4) . For radially long blades and variable root or tip 
radii, methods are proposed for analyzing the flow by a three- 
dimensional "throu^-flow" calculation, which considers the axial and 
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radial variations of the flow hut only a mean value in the circumferen- 
tial direction (references 5 and 6) . 0?he flow on the relative mean 
stream surface, which divides circumferentially the mass flow passing 
through the. channel between two blades into two equal parts, is taken 
to represent the mean flow through the blading. It is suggested in 
reference 5 that the shape of this mean stream surface and the correction 
factor b for a finite number of thick blades be obtained from the 
analysis of a number of two-dimensional flows along general surfaces of 
revolution starting at different inlet radii. A method was therefore 
developed at the MCA Lewis laboratory to determine the flow variation on 
these flow surfaces in a supersonic turbomachine. V?hen the flow surfaces 
may be approximated by cylindrical surfaces, the flow equations reduce -to 
the usual plane flow where the holograph characteristics are applicable 
for the flow analysis if the entrance shock is weak. 

TVi order to illustrate the effects of blade thickness and blade 
curvature and to determine, in general, the correlation between the 
shapes of the mean streamline apd the mean blade line , flow on the mean 
streamlines was determined for two supersonic cascades and a partly 
supersonic hyperbolic nozzle. 

This method, in addition to its use for flow analysis of a given 
blading, can be used to design blade elements along stream filaments of 
revolution having arbitrary thickness variation along the flow path. 


SYMBOLS 

The following symbols were used in this report; 

A area normal to velocity 

a local speed of sound 

b correction factor for finite number of thick blades 

c line tangent to characteristic curve 

Cp specific heat at constant pressure 

h static enthalpy 

I h + 

J,K,L coefficients of partial derivatives of ijr 
M,N 
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7,cp 

M 


m 


orthogonal coordinates on mean surface of revolution 
local Mach number 

mass flow between suction surface of blade and streamline, 

h 

pw dE 


P 

R 

r 

s 

T 

t 

W 

P 

7 

X 


P 

a 

T 


blade pitch or spacing 
gas constant 

radial distance from axis of machine (fig. l) 
entropy 

static temperature 

blade thickness in the y-dlrectlon 

resultant relative velocity 

rectangular coordinates for cylindrical flow with y-axls* 
chosen along line Joining leading edges of blade 

W 

-1 cp 

relative flow angle, tan ^ 
ratio of specific heats 

slope of characteristics on mean surface of revolution 


Mach angle, sin' 


■1 1 
M 


distance measured from suction to pressure surface normal to 
resultant velocity 

mass density 

slope of mean surface of revolution in meridional plane, 

normal thickness of stream filament of revolution 
stream function 
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cjo angular velocity of Hade 

Subscripts: 

i inlet condition 

meridional and circumferential components 
m on mean streamline 

r;Z components in r- and z-directions, respectively 

8 initial supersonic state 

t at throat of nozzle 

x^y components along x- and y-directions , respectively 

1,2 first and second family of characteristics, respectively 


METHOD 

On a general surface of revolution defined hy the orthogonal 
coordinates I and cp (figs. l(a) to 1(c)), the supersonic flow of a 
fluid along a stream filament of revolution of varying normal thickness 
T =» t(Z) is described by the-foHowlng forms of the flow equations. By 
analyzing the flow going in and out of the element shown in figure l(d) , 
the continuity equation for steady flow may be esqjressed by 


^(TpW^r) S(TpW^) 


”52 5 ^ 

When equation (l) is expanded and the relation 


0 


( 1 ) 




is used, equation (l) becomes 




= 0 


( 2 ) 
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For adiabatic flowj the entropy along a streamline remains constant for 
a nonviscous fluid j that is 


Ds .. Ss , ^ Ss ^ 


vhich allows the last term to he dropped from equation (2) . For a 
perfect gas. 


h => CpT 

and 


a^ 3 (y - l) h 


■When this substitution is made for a^ and the expression is multiplied 
1 

hy the eq'uation of continuity (2) becomes; 

Expressed in the coordinates I and cp, the equation of steady motion of 
a nonviscous fluid as given by equation (l4a) in reference 5 in the 
circumferential direction, is 


1 ^ 7 ^ 


ice + 

r 


2o3|sin a - ^ 


_2^ 


I Si 


T Ss 
r ^ 


r dcp 

From equation (3) , a stream function ijr is defined by: 


= 0 


H = ^TWxh 


1_ 

r-1 


(4) 


(5a) 


^ 


(5b) 


Substituting equations (5) into equation (4) results in the follow- 
ing relation: 
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2 


(irr + — 
r 


sin 


a + 


St , 


Si 




1 




sin a 



+ 



1-Sl _ T Ss \ 

r 5^ ” r 


1 Si I 1 St 
a2r ^ 


0 

( 6 ) 


If the symbols J, K; L, M, and IT are used to represent the coefficients 
of the partial derivatives in equation (6), it- may he rewritten as 


^ S^t , 2K S^t , L S^t , „ St . W St . . o 
•^^72 + TS ^-"; 2 ^ + + ° 


( 7 ) 


!The characteristics of equation (?) ar*e then 

J fj - (r f ) H. L . 0 

For equation (s), two real characteristics exist for supersonic 
flow, the slopes along which are 


( 8 ) 



K 

Vk^-jl 


l" 3rji = 

J “ 

J 

(9a) 

\ dZ/2 

K 

J 

(9h) 


When the polar coordinates transformation that p tan' 


Kn 

-1 JP 

W, 


and 


sin 


■"(h) 




sin' 




¥,2 + W(p2 


^6 used, the relation of X 


to the flow angle and the Mach angle is 
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= tan (p + p) 

(lOa) 


^2 “ "taJi (P - p) 

( 10 b) 

By the use of equations (?) and (9), the rate of change of 
i|f- derivatives along the characteristics with respect to Z is 
(see reference 5 for procedure) : Along c^^. 

obtained 


dZ \ SZ / r dZ \^J jSZ JrB^"^ 

( 11 a) 

and along cg 

d d ,N ^\|f ^ 

dz v^zj"^" dTl^/+ 

( 11 b) 


For ease in computation, it is more convenient to express equa- 
tions (ll) in terms of the magnitude of the resultant velocity W and 
the flow angle p as follows: Substituting equations (5) in equa- 

tions (ll) and differentiating yields the following expression along c-^: 


(^2 Wj 


- Wj +T 

a^ 


CpT o2 dZ 


^2 dZ " dZ ; 




9 W, - W) ^ + 
2 Z qj dZ 


1 {Ml - MW(p) + TWj sin a = 0 


Similarly, along 03 , 


^ ^1 ^Z ■ Hr)) — — + T f + (Xi W 7 - W ) ^ + 

^ ^ ^2 dZ ^ \ 1 dZ dZ / ^ dZ 


J (HWj - MW(p) + ^ TW^ sin C 7 = 0 


t 
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With 


and 


= W cos p 
Wjp = W sin p 


dZ " dZ 


2 2 j 


^ TT iW , 2 

dz - ^ dZ ^ 


and with tlie use of equations (lO), equations (ll) become 


1 w 
w 


tan ^ - tan^ 


03 r sin a 


a2 


— + i ^ - 

dz T dZ ^ cos p 


■2 sin a cos p 

sin 3) 


^ 1 N COB 3 - M Bln P 
J (Ag cos 3 - sin 3) 


(12a) 

(iaz3) 


03^r 


sin tj , 1 dl . 1 dT 


a“ 


+ -" i :tl + g cos 3 ^ ^ H cos 3 - M Bln 3 

^ 3T T dZ r(Aj_ cos 3 - sin 3) J cos 3 - sin 3J 


along Cj_ and C 2 , respectively. 

From two points a and b (fig. 2) where the flow is known^ the 
tangents to the. characteristics curves at these points are determined 
from equations (9a) and (9b); respectively; and intersect at a point c. 
By writing equation (l2a) (equation along cj) at a and equation (l2b) 
(equation along cg) at b in finite difference form and solving simul- 
taneously; the new flow values at c are detenninedj thuS; along C]_ 
at a; 


1 

Wa 



tan Pa, 



p 

tan Pa 


Oirx sin 
^2 


(13) 


and along at b; 


Wb Zc - Zb 


tan pi 3 


Pc - Pb 

^c ~ ^b 


o 

„ co“r sin CF^ 
tan^ p^ E + 

ay 


(14) 


The unknowns in both equations (l3) and (l4) are W^, and AZ is 

determined by the intersection of the tangents to the characteristics 
from points a and b so that - Zq and Z^ - Z-^, are known in the 
equation. 
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For a cliaracterlstic curve that Intersects the passage houndary^ 
such as is shown at d or e in figure 2, the situation is a little 
different. In a direct problem with a given blade configuration, the 
angle p at the intersecting boundary point is known. In the inverse 
or blade-design problem, either p at the Intersecting boundary point 
is known from a desirable turning specified at that point or the magni- 
tude of the resultant velocity W at the intersecting boundary point is 
known from a specified velocity distribution on the blade. In all cases, 
either equation (l3) or (l4) written for the characteristic from 
d or e determines the unknown W or p at the Intersecting boundary 
point, whereas, as shown previously, two simultaneous equations along 
intersecting characteristics are required inside the passage. In 
general, after W and p are determined at a new point, other fluid 
properties at that point are calculated and the characteristic directions 
constructed to repeat the process progressively downstream. (For irrota- 
tional inlet flow with uniform I and s and when the rotationality 
Introduced by the entrance shocks is neglected, I and s are uniform 
throughout for adiabatic flow, and the derivatives involved in equa- 
tions (l2) become zero. Otherwise, the derivatives are to be evaluated 
by using the values at the inlet on different streamlines, the constancy 
along the streamline before and after shock, and the changes across the 
shock. ) 

Thus a method is available either to analyze the blade- to-blade 
flow variation along a given arbitrary stream filament of revolution in 
a supersonic turbomachine, or to design the blade element on a given 
arbitrary stream filament of revolution for a specified turning distri- 
bution or a specified velocity distribution around the blade. The con- 
figuration of the stream filaments (fig. l(a)) is taken from a 
through-flow calculation (references 5 and 6) . For a direct problem, 
successive ceilculations between the through-flow calculation and the 
present calculation are necessary until the solution converges. The 
through-flow calculation gives the configuration of stream filaments of 
revolution, and the calculation on these stream filaments of revolution 
gives the shape of the mean stream sinrface and the factor b used in 
reference 5, which accounts for the effect of blade thickness and 
curvature. In an Inverse problem, the calculation is shortened if the 
estimated values of the factor b used in the through-flow calculation 
for a desirable blade-thickness distribution results in a blade dimension 
which has a thickness distribution close enough to the desirable one. 

In order to give a general idea of the order of magnitude of this 
factor b for the effect of the blade thickness and curvature on the 
mean flow and the closeness of the shape of the mean stream surface to 
the blade mean surface, an analysis is made for the simpler case of 
cylindrical flow utilizing available data for two two-dimensional 
high- solidity cascade recently • investigated by Liccinl at the MCA 
Langley laboratory (unpublished data) and for a partly supersonic nozzle 
Investigated by Emmons (reference 7) . 
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AMLYSIS OF FLOW IK TWO SUFERSOHIC CASCADES AKD 


OKE PARTLY SUPERSONIC NOZZLE 

The shapes of the two high- solidity 90° turning- angle cascades are 
shown in figure 3. Each caacade has an inlet Mach number of 1.78 and 
a 10° wedge angle at the leading edge. (The blade is so designed that — 
the leading-edge shock is cancelled after the first reflection.) The 
shape of the partly supersonic nozzle is shown in figure 4. (Because 
the nozzle is symmetrical, only half of it is shown.) 


In the case of the symmetrical nozzle, the center line is the mean 
streamline. For the two 90° cascades, the mean streamline (the one 
which divides the mass flow in the channel into two equal parts) has to 
be established from the data given on the characteristic net of the 
blade design. The mass flow m at any point in the channel is deter- 
mined by Integrating across the channel the local specific mass flow 

ri , 

given by the characteristic net} thus, m = j pW d&. The- values of 


o 

the local specific mass flow can be obtained, by using the table given in 
reference 8 or similar tables. The shape of the mean streamline is to 
be compared with that of. the mean blade line, which is obtained by 
taking the mean of the blade coordinates in the pitch direction. 


The value of the specific mass flow W or W^ is then obtained 
along the mean streamline downstream of its intersection with the 
leading-edge shock. The discontinuity in the flow ciirves at the- point 
where the mean streamline Intersects the "reflection" is smoothed out 
in order to emphasize the general trend. These values are to be com- 
pared with the corresponding variation in the channel width (or area) 
perpendicular tcr- W or W^ in order to determine if there is any corre- 
lation between the variations in specific mass flow on the mean stream- 
line and the corresponding area. 


In the case of. the nozzle, the flow variation on the mean stream- 
line is obtained from, the Mach number given in reference 7. The vari- 
ation in specific mass flow is coit^jared with the channel-area ratio. 

The results obtained in the two 90° turning passages are given in 
figures 3 and 5 to 9. The resulta obtained with the symmetrical, nozzle 
are shown in figures 4^- 10, and 11. 

In f-lgures 3(a) and 3(b), the mean streamline is compared with the 
mean channel line. In figure 3(a), where the cascade is of the higher 
solidity and has thinner blades, the shape of the mean streamline - 
closely approximated that of the mean camber line} however, in the 
second cascade (fig. 3(b)), which has a lower- solidity and thicker 
blades, the mean streamline deviates more from the mean blade line. 
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The specific mass flow on the mean streamline based on the resultant 
velocity as shown in figure 5 indicates clearly that the effect of blade 
thiclnaess and cTxrvature on the throiogh flow is to increase^ in both 
cases Investigated, the specific mass flow along the mean streamline a 
chordwlse average of 9 percent over that given by the area reduction 
(due to blade thickness) based on a one- dimensional calculation. If a 
comparison of specific mass flow in the axial direction with variation 
of channel width in the y-direction is made (fig. 6), as is done in 
reference 9 for subsonic flow, corresponding trends are obtained, and 
the increase in specific mass flow again averages about 9 percent 
chordwlse over the area ratio. This increase is more than two times the 
increase obtained in reference 9 . 

The variations of axial and tangential velocity components on the 
mean streamline are shown in figures 7 and 8, respectively. The vari- 
ations of the velocity across the channel in the two cascades at the 
stations indicated on figure 3 are shown in figure 9. 

The result obtained in the nozzle problem is shown in figures 10 
and 11. The results shown in figure 10 correspond to the nozzle shown 
in figure 15 of reference 7, which includes a small portion of super- 
sonic flow in the divergent portion of the nozzle; whereas figure 11 
corresponds to figure 16 of reference 7, with supersonic flow extending 
throughout the divergent portion and the Mach number increasing 
beyond 2. For the subsonic flow in the convergent portion, the local 
specific mass flow and area are noimalized by the values at the throat; 
for the flow in the divergent portion (partly supersonic, fig. 8, and 
purely supersonic, fig. 9), the local specific mass flow and area are 
normalized with respect to the conditions where the flow is Initially 
supersonic. The comparison between the specific mass flow and area 
ratio is similar to the previous result. The increase in the specific 
mass flow on the mean streamline is entirely due to the area reduction 
(no turning) and is about 8 percent higher than a one- dimensional 
correction. 

The results obtained in these analyses give some idea regarding the 
closeness between the shape of the blade mean line and mean streamline 
and of the order of magnitude of the blade thickness and curvature 
effect for through- flow calculations. For an approximate through- flow 
solution of a similar blading, these results may be directly applied. 

For different types of blading, similar analyses can be made for the 
two-dimensional flow on cylindrical surfaces or surfaces of revolution, 
which are formed by fluid particles initially lying on circular arcs at 
a number of radii upstream of the machine by using the general relations 
derived herein. The mean streamline obtained at the different radii can 
be joined together to form the mean stream surface. The ratio of 
(p^x^m these calculations can be taken as the 

thickness correction factor b in the relations between stream function 
and velocities (equation (ill) of reference 5) . 
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More knowledge of this kind for a mjmber of-typical bladings will 
be very useful in the design of these bladings. From this knowledge a nd 
a certain blade -thickness distribution, which is desired from the blade 
strength and other considerations, a good estimate of the factor b can 
be made and used in the through-flow calculation. From this solution, 
the blade section on a number of surfaces of revolution can be designed 
by the method given in the report. 


SUMMARY OF RESULTS 

Analysis of the supersonic flov in two two-dimensional 90° turning 
passages of high solidity and in a partly supersonic symmetrical nozzle 
shows that there is, in general, significant deviation of the mean 
streamline shape from that of the mean blade' line and 'that the effect of 
blade thickness and blade curvature on the specific mass flow along the 
mean streamline is to increase the specific mass flow along the mean 
streamline. about 9 percent above that given by a one- dimensional 
estimate. In order to determine these effects more accurately for 
turbomachines of arbitrary hub and casing shapes to be used for the 
three-dimensional through- flow calculations, a method is developed for 
the-determination of the supersonic flow along stream surfaces of revo- 
lution in turbomachines. In this method the shapes of the stream 
surfaces are arbitrary, and the method also takes into account the 
distance between adjacent- stream surfaces, which distance varies along 
the flow path. Thus, the method can be applied to turbomachines with 
arbitrary hub and casing shapes. 

In addition to their use for direct problems, these equations can 
be used to design blade elements in supersonic flow along an arbitrary 
stream filament of revolution in turbomachines. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 13, 1951 
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(a) StreamllrK in meridional plane . 



("b) Slream fUamen't of (®) Mean, stream tfnrfaoe 

revolution. revolution. 


Figure 1. - Flow on stream filament of revolution and 



(d) Element. 



face of revolution. 
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Figure 2. - Blade section and calculating points . 
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(a) Case a. 


Figure 3. 
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(a) Case a. (h) Case h. 

Figure 5. - Comparison between variations of specific mass flow on mean streamline 

and normal channel area ratio. 
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Chord, percent 

(a) Case a. (b) Case b. 


Figure 6, - Comparison of specific mass flow in x-dlrection with channel area ratio In 

y-direotlon. 
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Figure 7 . - Variation of axial veloolty on mean 

streamline . 
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Figure 8. - Variation of tangential velocity on 
mean streamline. 



Plgu.i'e 9= - Velocity yarlatlon acrosa chaimel 
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Figure 10. - Comparleon of variation In mean specific mass flow and area ratio 
with short portion of supersonic flow in divergent section. 
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Figure 11. - Comparison of variation in mean specific mass flow and area ratio 
of nozzle with supersonic flow ttiroughout divergent section. 
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